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ONGOING GALACTIC ACCRETION: 



SIMULATIONS AND OBSERVATIONS OF CONDENSED GAS IN HOT 
HALOS 
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ABSTRACT 

Ongoing accretion onto galactic disks has been recently theorized to progress via the unstable cooling 
of the baryonic halo into condensed clouds. These clouds have been identified as analogous to the 
High- Velocity Clouds (HVCs) observed in HI in our Galaxy. Here we compare the distribution of 
HVCs observed around our own Galaxy and extra-planar gas around the Andromeda galaxy to these 
possible HVC analogs in a simulation of galaxy formation that naturally generates these condensed 
clouds. We find a very good correspondence between these observations and the simulation, in terms 
of number, angular size, velocity distribution, overall flux and flux distribution of the clouds. We 
show that condensed cloud accretion only accounts for ~ 0.2M Q / year of the current overall Galactic 
accretion in the simulations. We also find that the simulated halo clouds accelerate and become more 
massive as they fall toward the disk. The parameter space of the simulated clouds is consistent with all 
of the observed HVC complexes that have distance constraints, except the Magellanic Stream which 
is known to have a different origin. We also find that nearly half of these simulated halo clouds would 
be indistinguishable from lower-velocity gas and that this effect is strongest further from the disk of 
the galaxy, thus indicating a possible missing population of HVCs. These results indicate that the 
majority of HVCs are consistent with being infalling, condensed clouds that are a remnant of Galaxy 
formation. 

Subject headings: ISM: kinematics and dynamics, ISM: clouds, Galaxy: halo, galaxies: formation 



1. INTRODUCTION 

Galaxies have long been thought to form from the 
cooling of s hock-heated primord i al material in galacti c 
halos (e.g. iWhite fc Reesl fl978L IWhite fc Frenkl fl99lh . 
These galaxy formation models, wherein all of the gas 
inside some time-dependent "cooling radius" monolithi- 
cally collapses, cannot generate the continuous accretion 
need ed to reproduce the Galaxy's star-formation histor 
(e.g. iRocha-Pinto et al.ll2000D . iMaller & Bullockl (|200 
(hereafter MB04) showed that the process of the forma- 
tion of galaxies could be driven by a multi-phase cooling 
process, in which not all of the gas within the cooling 
radi us simply collapses (s ee also iKaufmann et al.ll2007l 
and ISommer-Larsenll2006h . but instead the hot (~ 10 6 
K) shock-heated gas condenses out into warm (~ 10 4 K) 
clouds which fall inwards and collide, contributing to the 
formation of the galaxy. As time goes on, clouds con- 
dense at lower densities, corresponding to greater dis- 
tances from the galactic center, and 'rain' inwards to- 
wards the disk. This formulation, wherein accretion is 
more continuous, avoids the "overcooling problem" that 
plagued monolithic collapse models. MB04 identified 
these condensations as being analogous to observed High- 
Velocity Clouds (HVCs). This raised the exciting possi- 
bility that observations of HVCs can be used as a direct 
measure of the rate of our Galaxy's growth. We now 
further consider that we may be able to study galaxy 
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formation quantitatively by studying various HVC ob- 
servables and mapping them to the various parameters 
of galaxy formation simulations. Such a course of study 
might help disentangle the large scale processes (cool- 
ing, merging, outflows, etc.) that shape today's galaxy 
population. 

HV Cs have been known for over 40 years (|Muller et al.l 
Il963[ ). having been first observed in the 21-cm hyper- 
fine transition of neutral hydrogen. They are observed 
to be moving at hundreds of km s _1 , beyond what 
would be expected on the basis of Galactic rotation, 
with an overall sense of infall. Large m aps have been 
made of the Galactic sky in HI (e.g. iPutman et al.l 
l2002t iKalberla et al.l [2005), providing a good under- 
standing of the gross observational characteristics of the 
HI HVCs. HVCs have also been observed in other 
atom ic transitions (e.g 
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Putman et al.ll2003a. 


Tufte et al. 



2002) and there are tentative detections of dust as well 
(|Mivill e-Desch §nes et al.ll2005l ). These observations are 
much more sparse than the HI observations, and do not 
give us a sense of the distribution of the entire HVCs pop- 
ulation, though they provide s ome useful fiducial char- 
acter istics, such as metallicity (|van Woerden fc Wakkerl 



2004) and rough estimates of ionization fraction (|Tuftd 
2004h . Though the halo has been shown to h ave signif- 
icant high ions toward some HVC sightlines (|Fox et al.l 
2006), the majority of Ha observations show limited ion- 
ization fractions in HVCs around 20%. 

The origin and physical character of HVCs is still not 
known. The Magellanic Stream, a complex of HVCs 
trailing the Magellanic Clouds, is certainly due to the 
interaction of the Galaxy and the Magellanic Clouds 
and may have come from some combination of ram- 
pressure stripping by the gaseous halo and tidal ef- 
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fects ([Putman et al.ll2~003af h Some Intermediate- Velocity 
Clouds (IVCs) have been shown to be very nearby (~ 1 
kpc from the Galactic disk) and have near-solar metal- 
licities, and are therefore thought to be part of a Galac- 
tic fountain, a formation scenario in which gas is kicked 
out of the disk and rains back down. The rest of the 
population of HVCs are thought either to be satellite de- 
bris akin to the Magellanic Stream or part of the cooling 
formation pr ocess. In the p ast both Galactic fountain 
models (e.g. Ide Aville3l2000t ) and Local Group models, 
wherein HVCs are dark-matter domi nated clouds hun - 
dreds of kpcs from the Galaxy (e.g. iBlitz et aT1ll999f h 
have also been invoked to explain HVCs. Neither of 
these models, though, is currently thought to be a suc- 
cessful explanation of the bulk of HVCs. Many HVCs 
have been shown to have distances beyond a few kpc 
(e.g. iThom et al.ll2006UWakkerll2001l). and low metallic- 
ities (for a reviewjyan Woerden fc Wakkenl200l . which 
conflicts with the fountain model predicti ons of nearby, 
disk- li ke gas. Observations of Andromeda (|Thilker et al.l 
12004 IWestmeier et al.l l2005j ) have excluded massive 
HVCs at great distances from that galaxy, which are a 
key prediction of Local Group models. 

In an effort to understa nd the significance and scope 
of a Galactic gaseous halo. ISommer-Larsenl (|2006T ) (here- 
after S-L06) ran numerical simulations of galaxy forma- 
tion in the ACDM cosmology. These simulations showed 
that in Milky- Way-sized galaxies, significant baryonic 
mass resides in the hot halo. S-L06 also showed that in- 
falling condensations are a natural consequence of such 
halos, and that it is possib l e to c ompare these clouds to 
observed HVCs. iPutmanl (|2006l ) found the total mass 
in HVCs is consistent with the mass in clouds found by 
S-L06. 

In this paper we study the condensations found in S- 
L06 as analogous to HVCs and compare them to the 
HI HVC population in our Galaxy in terms of their 
number, flux, spatial and velocity distribution. Since 
we do not expect the resolution of these simulations to 
be high enough to resolve the masses of all individual 
HVCs, we instead compare coherent groups of to HVCs 
to groups of simulated halo clouds. We will also com- 
pare the simulated condensations to the extra-planar gas 
observed around the Andromeda galaxy, which provides 
the opportunity to compare the projected distance of the 
clouds, albeit with much lower resolution and sensitivity. 
We show that these condensations are indeed consistent 
with the overall HVC population and examine the sim- 
ulation population for physical characteristics that may 
inform our future study of HVCs and ongoing galaxy for- 
mation. 

The structure of the paper is as follows: In ij^lwe dis- 
cuss observations of HVCs around the Milky Way and 
and extragalactic HVCs around the Andromeda galaxy, 
in we discuss the details of the S-L06 simulations and 
the halo clouds found in the simulations, in §4] we de- 
scribe the process by which we self-consistently associate 
clouds in complexes, in |JS]we compare the observational 
data with the simulation, in we discuss the implica- 
tions of a cooling-formation origin for HVC observables 
and we conclude in fJ7] 



2.1. Selecting HVCs from Galactic Observations 

As we are interested in the global properties of HVCs 
rather than their minutiae, which are beyond the reso- 
lution of the simulation, we wish to use a full sky sur- 
vey with consistent resolution and nomenclature. For 
this reason we use the updated Wakker an d Van Woer- 
den catalog (jWakker &: van Woerdenlll991l . hereafter the 
WvWcat). This catalog has the velocities , fluxe s and 
positions of more than 600 clouds; IWakkerl (|2004j) uses 
this same sample, and a more detailed description can 
be found therein. 

The WvWcat includes all clouds that have historically 
been catalogued as HVCs, and therefore there is some 
inconsistency in their kinematic selection. For consis- 
tency we exclude those clouds that do not fit a deviation 
velocity criteria and a classical Local Standard of Rest 
(LSR) velocity criteria. The deviation velocity, Udcv, is a 
measure of how the velocity of a cloud deviates from a 
relatively simple model of the Galactic disk. The model 
includes solid-body rotation in the Galactic cent er, a flat 
rotat ion curve at larger R and a flared disk (see IWakkerl 
120041 for details). This is a useful criteria, as acceptable 
Galactic Standard of Rest (GSR) velocities vary strongly 
as a function of angle on the sky, either requiring very 
harsh cuts (|wgsr| > 220 km s^ 1 ) or including a large 
amount of gas plainly part of the disk in the quadrants I 
and IV. Local Standard of Rest (LSR) velocities (which 
is to say heliocentric radial velocities corrected for the 
Sun's peculiar motion as compared to the mean motion 
of stars in the stellar vicinity) have been used historically 
to define the difference between HVCs and other Galac- 
tic clouds. We use the criteria |«dcv| > 60 km s" 1 and 
I'^lsrI > 90 km s^ 1 , both to excerpt the disk and to fulfill 
the classical definition of HVCs. These criteria are ap- 
plied to observed HVCs and simulated halo clouds alike. 
The |i>lsr| > 90 km s _1 has a very limited effect upon 
the selected complexes once the |vdcv| > 60 km s _1 is 
applied and our results do not significantly depend upon 
whether the |wlsr| > 90 km s _1 is applied. 

In addition to the velocity selection described above, 
we exclude clouds which have origins that are not cap- 
tured by the simulation, namely the Magellanic Stream 
and the Leading Arm of the Magellanic System. The 
Magellanic Stream and Leading Arm orig inated in the 
Magellanic Clouds (e.g. iBriins et al.ll2005T ) and though 
the simulation could capture a similar dramatic disrup- 
tion event of an accreting satellite, we find no such event 
in the simulated halo clouds (see section 13. ip . After ex- 
cluding these clouds and applying the velocity criteria, 
we are left with 412 clouds. 

The majority of our remaining 412 clouds in the WvW- 
cat are grouped into "complexes", or a large group of 
clouds in the same region of the sky with similar ve- 
locities. These complexes contain the majority of the HI 
flux of HVCs and, given the resolution of the simulations, 
form a useful comparison set. Unfortunately, these com- 
plexes were defined in a historical and qualitative manner 
that is impossible to reproduce for the simulated dataset. 
We therefore apply a quantitative, consistent method for 
generating complexes in both the simulated and observed 
data sets as described in section [U 



2. OBSERVATIONS 



3 



2.2. Extragalactic HVCs around the Andromeda 
Galaxy 

In addition to Galactic observations of HVCs, extra- 
planar gas, seemingly analogous to HVC complexes, has 
been discovered in other galaxies. In particular, recent 
observations of Andromeda (M31) have s hown a large 
numb er of distinct extra-planar clouds (|Thilker et al.l 
l2004h (hereafter T04). To avoid confusion, we will refer 
to these clouds that reside in the halos of other galax- 
ies as extragalactic HVCs. The T04 observations are 
the most comprehensive study of extragalactic HVCs and 
have enough sensitivity to detect a large fraction of con- 
densed clouds analogous to those from the S-L06 simula- 
tions. These observations have an advantage over Galac- 
tic observations in that they have a projected galacto- 
centric radial distance for each cloud, as well as a rel- 
atively accurate distance and therefore HI mass for the 
clouds. At a distance of 775 kpc Andromeda is the near- 
est spiral galaxy to our own Milky Way. With a mas s 
comparable to that of the Milky Way (jSeigar et alj2 006) , 
and without evidence of recent major mergers, we expect 
Andromeda to have a relatively similar recent formation 
history, and therefore a similar population of halo clouds, 
to our own Galaxy. The T04 observations cover a 94 kpc 
x 94 kpc square at Andromeda which, though smaller 
than the simulation domain, overlaps with the bulk of 
the simulated clouds. The resolution is 2 kpc with ca- 
pacity to detect clouds down to a few xl0 5 M Q of HI, 
depending upon their size, which is comparable to the 
mass resolu tion of the simulation . These observations 
(along with IWestmeier et al.l l2005f ) show that there ex- 
ists a significant population of extragalactic HVCs within 
50 kpc of Andromeda's disk, with masses ranging from 
the sensitivity limit up to 10 7 Mq, and that there are 
not extragalactic HVCs with masses > 10 6 M© outside 
R = 50 kpc; these are significant constraints to which we 
can compare the simulations. 

3. SIMULATIONS 

The code used for the simulations is a significantly 
improved version of the TreeSPH code, which has 
been used previously for ga laxy formation simulations 
(|Sommer-Larsen et al.l [2003) . The main improvements 
over the previous version are: (1) The "conserva- 
tive" entropy equation so lving scheme suggested by 
ISpringel fc Hernqu ist ( 2002) has been adopted. (2) Non- 
instantaneous gas recycling and chemical evolution, trac- 
ing 10 elements (H, He, C, N, O, Mg, Si, S, C a and Fe), 
has bee n inc orporated in the code following iLia et af] 
(|2002bh and ILia et alJ (|2002ah : the algorithm includes 
supernovae of type II and type la, and mass loss from 
stars of all masses. (3) Atomic radiative cooling depend- 
ing both on the metal abundance of t he gas and on the 
meta- galactic UV field, modeled after lHaardt &: Madaul 
( 1996) is invoked, as well as simplified treatment of ra- 
diative transfer, switching off the UV field where the gas 
becomes optically thick to Lyman limit photons on scales 
o f ~ 1 kpc. 

ISommer-Larsenl (|2006h selected a Milky- Way-like 
galaxy from a cosmological simulation and simulated its 
gaseous halo at extremely high resolution. The purpose 
of the experiment was to establish how large a mass 
fraction of the hot gas halos, shown to contribute sig- 



nificantly to the baryonic mass budget of such galaxies, 
condense into "warm" (T~10 4 K) clouds by thermal in- 
stability. The result of the experiment was that only a 
few percent of the hot gas mass forms warm clouds; it 
was suggested that these clouds would be the equivalent 
of the HVCs. For the purpose of this paper, specifically 
addressing the properties of these warm clouds, an en- 
larged version of the above simulation was performed. 

The base of the experiment was a 3.2 xlO 5 particle, 
fully cosmological simulation of a disk galaxy, which 
at z=Q has a characteristic circular velocity of V c =224 
km/s, very similar to that of the Milky Way. At i=10.0 
Gyr (z~0.3), all gas particles within 250 kpc galacto- 
centric distance are split in eight particles of mass 1 / 8th 
the original value and gravity softening length (inverse 
gravity force resolution) 1/2 of the original value. The 
simulation is then continued for 200 Myr, and then all 
gas particles within 100 kpc galacto-centric distance are 
again split in eight particles of mass 1/8 and gravity soft- 
ening length 1/2 of the previous values. At this point the 
simulation totals l.lxlO 6 particles, of which 8.6xl0 5 are 
gas particles. The warm/hot gas in the inner 100 kpc of 
the galaxy halo is then resolved with particles of mass 
m g as=11700 Mq and gravity soften ing length 128 pc. As 
discussed by ISommer-Larsenl (|2006|) this enables the sim- 
ulation to resolve HVCs down to masses of ~ 3 x 10 5 Af Q 
within 100 kpc galacto-centric distance. The (by now) 
very high resolution simulation is run for 500 Myr, un- 
til i=10.7 Gyr, at which point it was terminated due to 
the heavy computational load. By i=10.3 Gyr more than 
100 condensed halo clouds have formed, with the number 
increasing slightly over the next 300 Myr. All HVCs are 
found at galacto-centric distances less than 100 kpc, no 
HVCs were found in the region 100<r<250 kpc down to 
the mass resolution of ~ 3 x 10 6 M Q in this region of the 
halo (Sommer-Larsen 2006). We note that the necessary 
condition for onset of thermal instability, viz. t\ < T coo \, 
is satisfied everywhere in the hot halo gas (t\ is the sound 
crossing time, which is taken to be ~ 2/isph/c s , where 
^sph is the local SPH smoothing length and c s is the 
sound speed; t coo i = is the timescale for radiative 
cooling). 

A number of physical processes can prevent the for- 
mation of warm (T ~ 10 4 K) clouds out of (T - 10 6 
K) gas by thermal instability, or lead to the destruc- 
tion of warm clouds moving in such hot gas. In partic- 
ular, as discussed by MB04, these include thermal con- 
duction, Kelvin-Helmhotz (KH) instability and conduc- 
tive evaporation. This will impose a lower limit on the 
warm cloud mass, below which the above processes sig- 
nificantly affect the clouds before they are accreted onto 
the disk of the galaxy. As the number of SPH parti- 
cles in the warm clouds are typically < 500 (section I3.1|) 
and the simulation was run with zero thermal conduc- 
tivity, it is important to demonstrate that the outcome 
of th e simulation is not a ffected by numerical resolu- 
tion (|Sommer-Larsenll2006f ) or lack of realistic modelling 
of the physics involved. Moreover, lAgertz et al.l (|2006f ) 
showed that the SPH method may have fundemental 
probl ems in mode l ling t he KH instability. On the other 
hand, IVietri et al.l (|1997| ) showed that KH cloud- breakup 
by the dominant "champagne effect" is highly suppressed 
by radiative cooling, which is included in the present sim- 
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ill ations, but not in the numerical experiments performed 
bv lAgertz etafl (|2006f ). 

As discussed in i j3.2l and H6.41 the total mass of the 
warm cloud system in the simulation is ~ 10 8 Mq, and 
the accretion rate onto the disk is ~ 0.2 M /yr. This 
gives a characteristic accretion timescale of ~ 5 x 10 8 yr. 
For such an accretion timescale, it follows from the work 
of MB 04 that warm clouds more massive than ~ 10 4 M 
will not be affected by thermal conduction/evaporation 
or KH breakup. As discussed in section H3.2I the masses 
of the warm clouds identified in the present simulation 
lie in the range ~ 10 5 — 5 x 10 6 M , so even at much 
higher numerical resolution and with thermal conduction 
included, warm clouds in the above mass range should 
neither be prevented from forming nor suffer significant 
destruction. 

Based on the evolution of the underlying, lower mass 
resolution simulation from t ~ to t=10 Gyr it appears 
that some of the halo clouds may have been seeded by 
the remains of "cold accretion" events taking place much 
earlier in the history of the galaxy, as discussed by SL06. 
However, given the much lower resolution of the main 
underlying galaxy simulation, starting at Zi=39 and run- 
ning to £=10.0 Gyr, it is not possible to give a detailed 
discussion of the origin of the halo clouds on the basis 
of the present simulations. These constraints also im- 
ply that the results obtained in this paper should, in 
general, be regarded as preliminary. Eventually, simula- 
tions starting at early times and of yet higher resolution 
should be undertaken, although such simulations are un- 
fortunately computationally prohibitive at present. 

3.1. Selecting individual HVCs from the simulation 

We identified potential halo HVCs in three snapshots 
at 300, 400 and 600 Myr (t=10.3, 10.4 and 10.6 Gyr). 
First, all "seed" SPH particles in the halo, satisfying 
tih > nH.trig and T<3x 10 4 K, were identified. Second, a 
gas particle group finder was used to identify all SPH par- 
ticles in coherent regions in the halo, surrounding these 
"seed" particles, and satisfying Uh > nH,min- Third, 
only SPH particles in these regions satisfying T<3xl0 4 
K were retained. It is found that with n#, tr i 9 ~10 -2 cm -3 
and n# im j n ~10 -3 ' 5 cm -3 one identifies neutral or partly 
photo-ionized T~l-3xl0 4 K gas in HVCs and satel- 
lite galaxies (SPH particles in the coherent regions of 
T>3x 10 4 K typically have T~10 6 K and are almost fully 
collisionally ionized) , hence these density thresholds were 
adopted. Subsequently, the 7 satellite galaxies identified 
around this galaxy are removed on the basis of these sys- 
tems containing (1) gas of high central density (n#>l-10 
cm -3 ), (2) stars and (3) dark matter. We note that none 
of these 7 satellite galaxies are accompanied by a stripped 
tail of gas analogous to the Magellanic Stream in any of 
the 3 snapshots. Indeed, there is no correlation between 
the location in position-velocity space of the satellites 
and the selected condensed clouds. We conclude from 
this that though it is possible that some of the condensed 
clouds we see were at one point associated with a small 
satellites, there is no evidence for halo clouds having been 
stripped from such satellites recently. We now describe 
the characteristics of these simulated halo clouds. 

3.2. Physical characteristics of simulated halo clouds 



The three snapshots in the simulation (300, 400 and 
600 Myrs) have 113, 128 and 130 identified halo clouds, 
respectively (hearafter, when describing a character- 
istic of the simulation, three quantities in a succes- 
sion will refer t o the se three snapshots respectively). 
ISommer-Larse n (2006) describes the mass distribution 
of the clouds in the simulation, which is consistent from 
snapshot to snapshot. Simulated halo cloud masses range 
from 10 5 Mq to 5 x 1O 6 M , with total masses ranging 
from 8.8 x 1O 7 M to 10.7 x 1O 7 M . Note that in the sim- 
ulation of SL06, the highest resolution was achieved only 
in the region out to 50 kpc galactocentric radius, as com- 
pared to 100 kpc for the present simulation. Compared 
to the former simulation the number of warm clouds in- 
creased by about 50%, but the total mass in warm clouds 
by less than 25%. Also note that this total mass is much 
lower than the - 2 x 1O 1O M proposed in MB04. Fig- 
ure [T] shows the distribution of halo clouds in the three 
snapshots, indicating both masses and velocity vectors. 
Note that the halo clouds in each snapshot are not inde- 
pendent; a typical cloud moving at 100 km s" 1 will only 
traverse 30 kpc from the first to last snapshot, so large 
structures that are consistent from snapshot to snap- 
shot may indeed be related. The cloud distributions are 
roughly spherically symmetric and have an overall sense 
of infall. Clouds tend to have velocity vectors similar to 
their neighbors and show large-scale (up to 30 kpc) in- 
flow structures. This is worthy of note as large, coherent 
structures in HVCs are sometimes cited as qualitative 
evidence of a satellite-accretion origin. Figure [2] shows 
the number density and mass density profiles of each of 
the snapshots. We note that the 300 Myr snapshot has a 
noticeably lower number and mass density of clouds near 
the disk. This effect is consistent with not enough time 
having elapsed since the beginning of the highest reso- 
lution simulation for the most massive clouds to have 
formed and fallen into the center of the potential. The 
400 Myr and 600 Myr snapshots have very similar den- 
sity structures, indicating a converged and continuous 
halo cloud lifecycle. 

In Figure [3] we plot the distribution of speed in the 
galaxy rest frame with respect to galactocentric radius 
of the halo clouds for each of the three snapshots. The 
most striking feature of this plot is the increasing velocity 
of clouds with decreasing radius, consistent over time, 
which is also visible in Figure [TJ This should not be a 
terribly surprising result, as clouds deeper in a galactic 
potential will typically have more kinetic energy, but the 
simple idea that more distant HVCs may have velocities 
that are less extreme than other HVCs has not been much 
addressed in the study of HVCs to date. Also evident 
in Figure [3] is that at a given distance, more massive 
clouds typically move faster. This effect consistent with 
the supposition that drag from the gaseous halo on halo 
clouds has a significant kinematic effect, and that more 
massive clouds typically have higher column densities. 

3.3. Observing simulated halo clouds 

To make comparisons between the simulated halo 
clouds to observed HVCs, it is necessary to "observe" 
the simulated clouds to determine the angular position, 
flux and radial velocity of each cloud. We do this by 
positioning a "sun" in the plane of the simulated galac- 
tic disk, 8.5 kpc from the galactic center and measuring 
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Fig. 1. — The distribution of simulated HVCs in the three snapshots, showing HVC mass (left) and speed (right). Each snapshot is 

projected onto the x-y and y-z planes, and shows a ring of radius 15 kpc centered in the galactic plane for scale and orientation. Symbol 

area scales with the mass of the clouds at left, and symbol length scales with the speed at right. 



each of those quantities for each cloud in each simulation 
snapshot. Fluxes are determined by assuming the gas 
in HVCs is 70% hydrogen by mass (consistent with Big 
Bang nucleosynthesis), and that the clouds are optically 
transparent to 21-cm radiation. This transparency as- 
sumption is reasonable as HVCs have a peak brightness 
temperature o f a few K and a spin temp erature upwards 
of 1000 K (see lKulkarni k Heiled (| 1988ft for a discussion 
of the details of HI radiative transfer) . There remains an 
unknown neutral fraction in the clouds, as the simulation 
does not carefully track the ionization of the halo clouds, 
which scales linearly with flux. Radial velocities of halo 
clouds can be measured including (LSR frame) or exclud- 
ing (GSR frame) the effect of galactic rotation. There is 
an ambiguity in the angular position of the "sun" around 
the galactic center when taking these "observations" . To 
overcome this, we measure each parameter at 8 solar po- 
sitions: R = 8.5 kpc, <p — nn/4:, n — 1,2, 3. ..8, in galac- 
tic cylindrical coordinates, applying the velocity cuts as 
described in section 12.11 for each point of observation. 
When reporting an "observed" value we report the mean 
of the 8 "observations" , with an error that is the stan- 
dard deviation of these 8 "observations" . We report this 



measurement for the 300 Myr, 400 Myr and 600 Myr 
snapshots sequentially. As an example, the number of 
simulated halo clouds that would be considered HVCs, 
given our adopted radial velocity cuts, are 56 ± 7, 67 ± 13 
and 74 ± 9. 

4. CREATING HVC COMPLEXES 

To compare simulated halo clouds and observed HVCs 
it behooves us to group them into complexes in a con- 
sistent manner. The tendency for HVCs to be near one 
another in angle-velocity space (i.e. complexes) should 
be easily captured by the simulation, even if the total 
number of smaller clouds within a complex may not. One 
difficulty in generating a consistent method for determin- 
ing membership of a cloud to a complex is the differing 
limitations in sensitivity of the simulations and the obser- 
vations. The WvWcat is essentially a flux- limited catalog 
- very low-mass clouds that are particularly nearby are 
as easily observed as distant, massive clouds. The sim- 
ulation, by contrast, is limited in mass resolution. The 
sensitivity of the observations to resolve nearby, very low 
mass clouds may be the reason for the very large num- 
ber of low-flux clouds in the WvWcat that do not have 
analogs in the simulation data. To eliminate as many 
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Fig. 2. — The number density (top) and mass density (bottom) 
profiles for the simulated HVCs with galactocentric radius, in 10 
kpc bins. Each of the snapshots is plotted. There are no clouds 
between 90 and 100 kpc in the 600 Myr snapshot, as compared to 
1 and 3 clouds in the 300 and 400 Myr snapshots respectively. 



of these observed 'chaff' clouds as possible, we wish to 
determine the minimum flux an HVC can have in the 
WvWcat to have a mass and distance that can be cap- 
tured by the simulation. For clouds with fluxes of up to 
300 Jy km s^ 1 to have masses that can be resolved by 
the simulation, they must be at distances of at least ~ 50 
kpc. If these clouds did largely exist at these types of dis- 
tances and beyond, they would contribute > 4 x 1O 7 M 
of n eutral hydrogen t o hal o. The observations in T04 
and IWestmeier et al.l ()2005[ ) show that most of the mass 
of extragalactic HVCs (~ (3 — 4) x 10 7 M Q ) exists within 
50 kpc. Under the assumption that the Milky Way and 
Andromeda do not have qualitatively different HVC halo 
populations, clouds with fluxes below 300 Jy km s _1 can- 
not consistently be at distances greater than 50 kpc and 
therefore cannot be reproducible by the simulation. This 
assertion is consistent with H a observations that show 
many small clouds to have ionization indicating a dis - 
tance less than 40 kpc ()Putman et al.l2003bllTuftell2004D . 
We therefore make a flux cut at 300 Jy km s _1 , excluding 
all lower flux clouds in the WvWcat from our compar- 
ison. Though this cut removes 77 % of the individual 
HVCs in the WvWcat of interest, it removes only 5% of 
the flux. 

Now that we have reasonably comparable data sets, to 
determine which clouds should be associated into which 
complexes we define a 'distance measure' in position- 
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Fig. 3. — Three plots of HVC galactocentric radius versus speed 
(with respect to the center of mass of the Galaxy) of the three snap- 
shots in the S-L06 simulation. The size of the symbol scales as the 
square-root of the mass of the cloud. Clear trends of velocity with 
radius, mass with radius and velocity with mass at a fixed radius 
are evident. Six observed clouds and com plexes that have distance 
constr a ints are over-plotte d (d ata from [va n Wocrdcn & Wakkcr 
(2004), IWeiner et al.l <200j ) and IThom et* a l. (2006)), with arrows 
indicating where values are only upper or lower limits in distance 
or lower limits in velocity. 

velocity space between clouds: 



d = ve 2 + / 2 (Svy 



(i) 



Here, is the angular distance between two clouds, Sv 
is the difference in GSR velocity between two clouds and 
/ is a conversion factor; / parameterizes the significance 
we ascribe to the angle subtended by two clouds versus 
their difference in velocity in determining whether they 
are members of the same complex. We choose / = 0.5°/ 
km s _1 , broadly consistent with the clustering in histor- 
ical HVC complexes. We place two clouds in the same 
complex if their D is less than some number, Dq. A 
cloud with no neighbors is considered a complex of its 
own. The overall complex position and velocity are de- 
termined by a flux-weighted average of the constituent 
clouds. 

This formulation leaves us with a free parameter, Dq- 
As Dq is increased the number of complexes decreases 
roughly linearly for reasonable values of Dq in the simu- 
lated and WvWcat data sets, thus the data do not offer 
a specific scale at which cloud clustering takes place. We 
wish to choose Dq such that we maximize the identifi- 
cation of true complexes and minimize the identification 
of complexes that are the result of coincidental cloud su- 
perposition. To that end we scramble all of the positions 
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(1 and b) and velocities (fobs) of the simulated clouds 
within each snapshot, as well as in the WvWcat data, 
such that all of the coherent angular and velocity struc- 
tures are lost, while maintaining the distribution of each 
of these parameters individually. We then run the clus- 
tering algorithm on the scrambled and unscrambled data 
sets for all values of D . We assume that any complexes 
generated in the scrambled data sets are spurious and 
we find the average Do at which the difference in clus- 
tering, as parameterized by the number of complexes, is 
greatest between the true and scrambled data sets. This 
maximum occurs at Do ~ 25°, for each of the snap- 
shots as well as for the WvWcat data. We adopt this 
value for generating cloud complexes hereafter, and refer 
to complexes from the S-L06 halo clouds as 'simulated 
complexes' and complexes from the WvWcat HVCs as 
'observed complexes'. 

5. COMPARISON BETWEEN SIMULATED AND 
OBSERVED CLOUDS 

5.1. Milky Way HVC Complexes Compared to 
Simulated Cloud Complexes 

First we wish to compare the number and angular size 
of observed HVC complexes to the simulated complexes 
that fit our velocity selection criteria as outlined in [J2j 
Figure 2] shows each of the snapshots as observed from 
the Sun along with the observed HVC data. The number 
of complexes is similar, 20.7±1.5, 26.1±3.3and 29.0±2.0, 
as compared to 29 in the WvWcat. Some subtle evolu- 
tion exists in the number of clouds and complexes in 
the simulation data, consistent with the halo requiring 
more than 300 Myrs to reach an equilibrium in the con- 
densation of HVCs after the simulation has been run at 
high-resolution. Note that the simulated complexes have 
rather a random distribution on the sky, similar to ob- 
served complexes. Simulated complexes also range in 
angular size up to about 60°, consistent with range of 
angular sizes in the observed complexes. In Figure [5] 
wc show simulated complexes from four vantage points 
at the solar circle for each snapshot, demonstrating the 
variation in the distribution of complexes as a function 
of solar position. This overall agreement between simu- 
lation and observ ation is consistent with earlier results 

(pool . 



in IConnors et al 

We also wish to compare the velocities and fluxes of the 
simulated complexes and observed HVC complexes. Fig- 
ure [6] shows the fluxes and velocities for the simulated 
complexes in each of the snapshots from four vantage 
points around the solar circle, along with the observed 
HVC complex data. The total maximum HI flux in the 
simulated complexes are 4.5±1.9, 10.0±4.6 and 8.3±2.9 
xlO 5 Jy km/s, as compared to 5.7 xlO 5 Jy km/s in the 
compared sample of observed HVC complexes. Note that 
the maximum flux comes from the assumption that all 
of the simulated cloud mass is neutral, therefore if only 
a small fraction (~ 30%) of the simulated halo clouds 
were neutral, simulation and observation would have sig- 
nificant discrepancy in HI flux. The range in total flux 
within a single snapshot is due to the fact that ~ 50% 
of the overall flux comes from just a few simulated com- 
plexes closer than 10 kpc; as the observation point is 
rotated around the Galaxy, the distance to these com- 
plexes changes, changing the overall flux. This depen- 
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Fig. 4. — Four plots of HVC complex distribution on the Galactic 
sky. The first three are from the S-L06 simulations and the final 
plot is of the observed HVC complexes. Filled circles represent 
HVCs with the linear size of the circle scaling with the cube-root 
of the HI flux, as HI flux covers more than 3 orders of magnitude. 
In the simulation plots, clouds are gray if they do not fulfill the 
velocity selection criteria: [Udevl > 60 km s — 1 and I^lsrI > 90 
km s _1 . Clouds are colored by their velocities, and we assume 
an neutral fraction of 70% in the simulation plots. Lines connect 
clouds to the flux-weighted centers of the complexes to which they 
belong. Complexes are indicated by an empty circle, with fluxes 
scaling the same way with circle size. 

dence upon only a few large, local clouds, though consis- 
tent with observations, diminishes the usefulness of flux 
as a measure of the accuracy of the simulations. The av- 
erage velocities (GSR) are also similar: —81 ± 7, —82 ± 6 
and — 65±24 km s^ 1 , as compared to —85 km s^ 1 in the 
observed HVC complex dataset. The standard deviation 
of the distribution of simulated cloud complex velocities 
is 56 ±4, 68 ± 10 and 74 ± 25 km s _1 , as compared to 
75 km/s in the observed complex data set. These distri- 
butions of observables are consistent with the hypothesis 
that the condensed halo clouds in the S-L06 simulations 
are analogous to the Milky Way's population of HVCs. 

5.2. Andromeda's HVCs Compared with Simulated 
Cloud Complexes 

Here we compare observed quantities in the T04 An- 
dromeda dataset to our simulated data set. We make an 
analogy here between the observed extragalactic HVCs 
around Andromeda and simulated complexes, as the 
masses are similar. T04 find that there is ~ (3 — 4) x 
10 7 Mq in HI mass around Andromeda; we find 6.1, 
7.5 and 6.8 xl0 7 M© in our 3 simulation snapshots of 
300, 400 and 600 Myrs. This implies an neutral frac- 
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Fig. 5. — Twelve plots of Galactic 1 versus Galactic b for simulated HVC complexes. Columns from left to right are from each sequential 
snapshot (0.3, 0.4 and 0.6 Gyrs into the high-resolution simulation), and rows are from positions rotated by 90° around the solar circle. 
The symbol area scales with flux. Dark symbols represent simulated complexes and light symbols represent observed Galactic complexes, 
for comparison. Symbols are scaled using 70% neutral fraction. Note that the distribution of simulated complexes varies significantly with 
position of the Sun within the solar circle, but that the simulated complex distributions are all qualitatively similar to the observed complex 
distribution. 



tion ranging from 40% to 65% for the simulations to be 
consistent with the observations. The standard devia- 
tion of the projected velocities in the T04 Andromeda 
dataset is 126 km/s, although limiting the sample to ob- 
jects with masses greater than 5 x 10 5 M Q , where we ex- 
pect the sample to be complete, reduces this to 96 km/s. 
Also note that this velocity dispersion only pertains to 
the qualitatively-defined "objects" in the T04 dataset, 
which excludes diffuse gas closer to the Andromeda disk. 
The simulation snapshots, once projected to mimic the 
Andromeda viewing angle, have velocity standard devi- 
ations of 63, 71 and 89 km s _1 . 

The spatial distribution of the halo clouds is one of 
the most interesting quantities to pursue in this compari- 
son, as the Andromeda data set yields a minimum galac- 
tocentric distance for all observed extragalactic HVCs, 
a quantity not afforded by the Galactic data set. It 
is, unfortunately, a rather difficult quantity to extract, 
as the scale of distribution of the halo clouds is very 
sensitive to the sensitivity of the observations to small 
clouds, which is difficult to accurately characterize in 
these observations. We can conclude that both the An- 
dromeda data set and the simulation snapshots have de- 
tectable HVCs out to projected galactocentric radii of 



50 k pc, and lack clouds wi th M > 1O 6 M outside 50 
kpc (|Westmeier et al.l 12005). A more detailed analysis 
of the scale-length and HVC density profile, including 
lower mass clouds, will have to wait for deeper, higher 
resolution observations. 

6. ANALYSIS & DISCUSSION 

In this section we discuss the significance of ascrib- 
ing a cooling-formation origin consistent with the S-L06 
simulations to the observed HVC population for a few 
different HVC observables. We also discuss a few ways 
to observationally disprove the hypothesis that the bulk 
of observed HVCs have this origin. 

6.1. The Speed- Radius Domain 

Consistency between simulations and observations can 
be further examined through investigation of the speed- 
radius relationship for both populations. The distribu- 
tion that the simulated clouds present (see Figure [3|) has 
a very clear variation of Galactic rest-frame speed with 
galactocentric radius. The vast majority of the observed 
HVCs with distance constraints are consistent with the 
simulated HVCs. Only one observed complex in the 
Galactic sky with a known distance significantly strays 
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Fig. 6. — Twelve plots of GSR velocity versus HI flux for simu- 
lated HVC complexes and observed HVC complexes. In each plot 
the observed HVC complexes are depicted by gray filled circles for 
comparison. The dashes represent simulated complexes, with the 
right side of the dash representing 100% neutral and the left side of 
the dash representing 50% neutral. From left to right the plots are 
of 300, 400 and 600 Myrs into the simulation. From top to bottom 
the plots represent viewpoints of 0°, 90°, 180° and 270° around 
the solar circle. 

from this region: the Magellanic Stream. This discrep- 
ancy, of course, does not contradict the hypothesis that 
the rest of the HVCs are generated by a mechanism sim- 
ilar to the one that generated the S-L06 halo clouds; the 
Magellanic Stream was excluded from the analysis for the 
specific reason that is has a well-known, non-cooling ori- 
gin. The deviation of the Magellanic Stream from this 
area does indicate, however, that HVCs that are gen- 
erated from satellite accretion may occupy a different 
part of the radius-velocity parameter space and that it 
therefore may be possible to distinguish the origin of an 
individual cloud if these values can be measured. 

It is worth noting that the very high-velocity cloud 
(VHVC) 160.7-44.8-333 (annotated P07 in Figure EJ) has 
been suggested to h ave a galactocent ric distance of ~ 20 
kpc (jWeiner et al J 120011 : iPeek et al.ll2007t ). At a mini- 
mum velocity of 300 km s^ 1 GSR, it is also marginally 
outside of the domain of clouds that are generated in the 
S-L06 simulations. Further distance limits on VHVCs 
may help to show whether these clouds are broadly incon- 
sistent with cooling formation scenarios. Observations 
of VHVCs are particularly important as they have very 
large minimum Galactic rest frame speeds, and are there- 
fore plausibly excludable from the cooling-formation- 
origin domain. 



6.2. HVCs at Low Velocity 

A crucial byproduct of the analysis of simulated con- 
densed clouds as observable HVCs is that nearly half of 
all simulated clouds would not satisfy the observational 
velocity criteria for HVCs. Instead, these clouds would 
be construed as lower-velocity (and therefore nearby) 
gas. This is to say that if the bulk of HVCs are gen- 
erated by processes similar to those in the S-L06 simula- 
tions we are 'missing' as many HVCs as we observe and 
thus underestimate the total scale of HVCs by a factor 
of two. This effect increases with greater radius from the 
galaxy - the average velocity of clouds decreases with ra- 
dius (see Figure [3]), such that more distant clouds will be 
more often confused with local gas. The clouds with low 
radial velocity are startlingly asymmetrical. The ratio 
of the number of condensed clouds above \b\ of 50° with 
—50 km s _1 < Vqsr < km s _1 to clouds with km 
s _1 < V GS r < 50 km s" 1 is 7±4, 9±3 and 11±3 for each 
of the snapshots. If these clouds could be observationally 
deciphered from nearby LVCs and IVCs, this asymmetry 
could be a powerful test of this model for HVC forma- 
tion. It may be possible in the future to disentangle such 
clouds from local gas via their morphological character- 
istics, absorption to stars, low metallicity and very low 
dust-to-gas ratio. 

6.3. CHVC hypotheses 

Compact high- velocity clouds (CHVCs), are HVCs 
that are smaller than about 2° in projected size and are 
not directly linked with other HVCs in the sky. CHVCs 
have received significa nt attention in recent years (e.g. 
iBraun &: BurtorJ[l9 99) under the assumption that some 
or all of them are analogous to the observed large com- 
plexes of HVCs but at much greater distances. Were this 
true, CHVCs could dominate the mass of Galactic HVCs, 
and would be very important to understanding the struc- 
ture of the HVCs and the Galactic halo. This assump- 
tion hinges on the assumption that the mass of an HVC 
complex is relatively independent of distance from the 
Galaxy. If HVCs have a mass-distance distribution simi- 
lar to the simulated distribution, this assumption would 
prove false. The most massive halo clouds in the S-L06 
simulations are typically closest to the disk, and they 
do not have analogs 50 to 100 kpc from the galaxy. In 
addition to this, the simulated velocity distribution of 
local halo clouds (< 30 kpc) is not analogous to that of 
more distant halo clouds; distant halo clouds move slower 
in the S-L06 simulations. Thus if some CHVCs are in- 
deed distant objects (and the cooling-formation predic- 
tions are correct) there will exist a trend toward lower 
velocities in CHVCs. CHVC catalogs do n ot show a lower 
veloc ity dispersion than catalogs of HVCs (jPutman et al.l 
2002), but such a discrepancy may be hard to detect 
as small, local clouds could easily masquerade as dis- 
tant CHVCs (thus contaminating the sample) and disk 
gas may obscure a large fraction of the population of 
slow-moving, distant CHVCs. Conversely, it is possible 
to limit the contribution to the HVC population from 
cooling-formation HVCs if CHVCs with Vqsr > 200 km 
s _1 can be shown to be distant (R > 50 kpc). This is 
to say the velocity distribution of distant clouds may be 
a strong test of this cooling-formation scenario and the 
SL06 simulations. 
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Fig. 7. — The accretion rate onto the simulated galaxy 
as a function of radius, as defined by the equation M (R) = 
y-^i— n(R) M, r, ) ^ rp^g pj o ^_ s j lows ^- ne av erage across all three 
snapshots. Inside the dashed line the parameterization of the accre- 
tion rate with radius is no longer meaningful as clouds are merging 
with the disk. 

6.4. The HVC Accretion Rate 

In this section we compare the accretion rate in the S- 
L06 simulations, and the source of that accretion, to the 
accretion necessary to fuel the ongoing star formation ob- 
served in our Galaxy. In the S-L06 simulations, material 
cools out of the halo at different radii, allowing us to pa- 
rameterize the halo clouds accretion rate onto the galaxy 
as a function of radius in the simulations. At R=100 kpc 
the accretion rate is zero and the rate increases with de- 
creasing R. We define the halo clouds accretion rate at a 
given galacto-centric radius by 

i=n(R) y . . 

M(R)= E (2) 

1=1 

where Mj is the mass of a cloud, V{ is a cloud veloc- 
ity and fi is the radial unit vector at the cloud, n (R) 
is the number of clouds that exist in a spherical shell 
from R - dR/2 to R + dR/2. M (R) is not a well-defined 
metric within 10 kpc of the galactic center, as some of 
the halo clouds will be removed from the system when 
they collide with the disk, so we exclude HVCs within 
10 kpc from this analysis. We find that this halo clouds 
accretion rate, averaged over all three snapshots, can be 
easily fit by a line with M (0) = 0.22 ± 0.014M Q /year 
and dtil (R) /dR = -2.5 x 10~ 3 ± 2 x 10" 4 M Q /year/kpc 
(see Figure[7]). This disk accretion rate of 0.22 M Q /year 
is consistent with the assertion that HVC accretion is not 
the sole source of fuel for Galactic star formation unless 
the accretion r ate was higher at e arlier in the lifetime o f 
the Galaxy (see lPutmanll2006l and lSommer-Lar scn 2006). 
The monotonic and linear increase in accretion rate with 
decreasing radius demonstrates that the cooling process 
that drives halo cloud formation and growth operates at 
all radii within the halo in the S-L06 simulations. This 
result may indicate that HVCs form and grow through- 
out our own Galactic halo as well. This is contrary to a 
simpler picture wherein all HVCs form at some typical 
distance in the halo and fall toward the Galactic disk. 



7. CONCLUSION 

The number of complexes as identified in our TreeSPH 
simulation and their angular distribution on the sky is 
consistent with the Milky Way HVCs, when the veloc- 
ity selection effects and clustering effects are taken into 
account. The radial distribution of simulated condensed 
halo clouds is consistent with the T04 sample of extra- 
galactic HVCs around Andromeda and, in particular, the 
lack of massive (~ 10 6 M Q ) clouds at large radii (> 40 
kpc) in the simulation is consistent with observations of 
Andromeda which show a dearth of such clouds far from 
the disk. The flux and velocity distributions of simu- 
lated complexes are also consistent with observations of 
HVC complexes in the Milky Way and of extragalac- 
tic HVCs in Andromeda, and point to neutral fractions 
above ~ 30%. This lower limit is inconsistent with a 
"tip-of-the-iceberg" picture of HVCs, wherein HVCs are 
dominantly an ionized phenomenon and the observed HI 
is just a small fraction of the baryonic mass of the clouds. 
The physical and observational population characteris- 
tics of the simulated halo clouds are consistent over 300 
Myrs (excepting the nearest, most massive HVCs in the 
first snapshot) and are broadly independent of the point 
of observation chosen on the solar circle. We have found 
that simulated halo material condenses into clouds from 
R = 10 kpc to R =100 kpc in the halo, and that these 
clouds have an overall accretion rate of ~ 0.2M Q / year. 

We have also shown that the speed-distance domain 
is populated by simulated halo clouds only in a specific 
region, and that all HVCs with known distances and un- 
known origins (i.e. non- Magellanic) reside in this region. 
This points to a possible method for discriminating HVC 
origin given HVC speed. The simulated clouds are not 
always identifiable as HVCs, as they may have low pro- 
jected velocity. Indeed, we find that half of the simulated 
halo clouds do not show high enough projected velocities 
to be considered HVCs at all. Specifically, distant sim- 
ulated halo clouds are typically less massive and slower 
moving than nearby simulated halo clouds which makes 
them harder to distinguish from local gas. This implies 
that the populations of distant CHVCs and nearby HVCs 
are physically different, if we assume that CHVCs and 
HVCs are formed through a similar cooling mechanism. 

Observations using the Arecibo L-Band Feed Array 
will refine our understanding of both the Galactic HVC 
population with th e GALFA-HI observin g prog ram (e.g. 
IStanimirovic et alj 120061 and iPeek et alj|2007h and the 
distribution of HVCs around other galaxies w ith the 
AGES observing program (e.g. lAuld et al.|[2006f ). AGES 
will generate maps with similar fidelity to those in T04, 
but for more galaxies out to ~ 5 Mpc. Instruments com- 
ing online now, such as the Allen Telescope Array, will 
allow us to map the HI HVCs in the vicinity of other 
nearby galaxies with unprecedented efficiency and fur- 
ther compare these cooling-formation simulations to ob- 
served systems. Indeed, the full ATA-350 will be able 
to make maps comparable to those in T04 out to ~ 10 
Mpc in a single day's observation. In this way we may 
be able to determine whether there are characteristics of 
extragalactic HVC systems that can inform our under- 
standing of the formation of these galaxies and the varia- 
tion in the character of their baryonic halos. In the more 
distant future, the Square Kilometer Array will allow 
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us to extend this analysis to non-zero redshift, probing 
the cooling history of galaxies toward the age of merg- 
ers. To match these observational efforts, more detailed 
and comprehensive simulations will be required. To this 
end, cosmological simulations of even higher resolution 
than the ones described in the paper, and with a more 
sophisticated treatment of UVB and stellar ionizing pho- 
ton r adiative transfer (e.g. iRazoumov fc Sommer-Larsenl 
2007J) are underway. 
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